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<T^ One of the main goals of the COMPASS experiment at CERN is the 

i— h determination of the gluon contribution to the nucleon spin. To achieve 

this goal COMPASS uses a naturally polarised muon beam with an energy 
i — i of 160 GeV and a fixed polarised target. Two types of materials are used 

^ for the latter: 6 LiD (polarised deuterons) during the years of 2002-2006 

^ and NH 3 (polarised protons) in 2007. The gluons in the nucleon can be 

dj accessed directly via the Photon Gluon Fusion (PGF) process. Among 

r^, the channels studied by COMPASS, the production of charmed mesons is 

the one that tags a PGF interaction in the most clean and efficient way. 

This talk presents a result for the gluon polarisation, Ag/g, which is based 

on a measurement of the spin asymmetry resulting from the production of 
^ D° mesons. These mesons are reconstructed through the invariant mass of 

their decay products. The statistical significance of the D° mass spectra 
OO has been improved significantly using a new method based on Neural 

Networks. The Ag/g result is also presented using a next-to-leading order 
y—{ (NLO-QCD) analysis of the jiN — > qq reaction. Such correction is relevant 

and was for the first time applied to an experimental measurement of the 

gluon polarisation. 
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1 Introduction 



It is known since the twenties that the nucleons carry half integer spin. However, 
what still remains unknown to the present day is the distribution of the 1/2 spin 
among the nucleon constituents. Until the eighties it was believed that the nucleon 
spin was distributed among its three valence quarks. In fact, assuming an SU(2) 
flavour symmetry and an SU(2) spin symmetry we obtain the following expression for 
the static wave function of the proton (the spin projection is considered to be parallel 
to the quantisation axis): 



|p|> = "7= t 2 \ututdi) - \u^uld^) - |«4,utdt) + (u d) ]. (1) 
v 18 

Eq. (1) allows us to calculate the contribution of the u and d quarks to the proton 
spin: 

Au=(pt\N u t -N u i |pt) = ^(10 - 2) = ^. (2) 

Ad=(pt\N d t -7Vabt> = ^(2 - 4) = -^. (3) 

where N q f (N q I) represents the number of quarks of flavour q with a spin projection 
parallel (anti-parallel) to the proton spin. Therefore, we obtain for the total helicity 
of quarks: 

AS:=A M + Ad = (| - |) = 1- (4) 

Indeed, from Eq. (4) we conclude that all the nucleon spin is carried by the va- 
lence quarks. However, in 1988 the European Muon Collaboration (EMC) at CERN 
astonished the scientific community by publishing a surprising result for the total 
contribution of quarks to the nucleon spin (AS). The experimental value obtained 
by EMC was very small, and even compatible with zero [1]: AS = 0.12 ± 0.17. This 
result gave rise to the so-called spin crisis of the nucleon because it could not be 
reconciled with the theoretical predictions: even applying the necessary relativistic 
corrections [2], the expectation value for AS is about 0.60 [3]. Since then, polarised 
lepton-nucleon scattering experiments were performed at CERN by SMC [4] and 
COMPASS [5], at SLAC [6], at DESY [7] and at JLAB [8] as well as in polarised 
proton-proton collisions at RICH [9, 10]. The goal is to extract the parton helicity 
distributions in the nucleon, using a perturbative QCD analysis. The contribution of 
quarks to the nucleon spin is nowadays confirmed to be 30% (AS = 0.30 ±0.01 ±0.01 
[5]). The fact that this number is still quite below the expected one leads us to the 
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following question: where is the remaining part of the nucleon spin? In QCD the 
nucleon spin projection (in units of H) may be decomposed into the quark and gluon 
helicities, AS and AG, and also into their orbital angular momenta L q and L g : 

\ = J q + J g = QAE + L 9 ) + (AG + . (5) 

Knowing that the gluons were the solution to the problem of the missing momentum 
in the nucleon, the obvious approach to solve this spin puzzle (Eq. (5)) would be 
the determination of the gluon helicity AG (note that AG = Jq Ag(x)dx). This was 
a strong motivation for measuring the gluon helicity distribution as a function of the 
fraction x of the gluon momentum, Ag(x), in dedicated experiments like COMPASS. 
Another motivation is found on the fact that in QCD the U(l) anomaly generates a 
contribution from gluons to the measured singlet axial coupling a (Q 2 ). This implies 
that AS is scheme dependent and may differ from the observable ao(Q 2 ), while AG 
is scheme independent at least up to the NLO approximation. In the AB (Adler- 
Bardeen) scheme, AS is independent of Q 2 [11]. It is related to the MS (Minimum 
Subtraction) scheme through the following equation: 

AS^(Q 2 ) = AS AB - n f as j^AG(Q 2 ). (6) 

Assuming the Ellis- Jaffe prediction of AT. AB 0.6, very large values of AG(Q 2 ) are 
required in order to reconcile the theoretical prediction with the experimental result: 
AS^(Q 2 ) = a (Q 2 ) » 0.3. 

Due to a limited range in Q 2 covered by the polarised experiments, the QCD evolu- 
tion of the DGLAP equations shows a very limited sensitivity to the gluon helicity 
distribution, Ag(x), and to its first moment, AG. Therefore, the determination of 
Ag{x) from QCD evolution must be complemented by direct measurements. The 
average gluon polarisation in a limited range of x, (AG/G) X , has been determined in 
a model independent way, from the photon gluon fusion process, by HERMES [12], 
SMC [13] and COMPASS [14, 15]. In this work we present the latest COMPASS 
results obtained from an open charm analysis. Basically, the PGF process is tagged 
by the detection of the decay products of charmed mesons. 



2 The COMPASS Experiment 

COMPASS is a fixed target experiment located in the CERN North Area, at the 
end of the M2 beam line of the Super Proton Synchrotron (SPS) accelerator. In 
order to study the gluon polarisation in the nucleon, COMPASS uses a naturally 
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polarised muon beam and a fixed polarised target. The average beam polarisation 
in the laboratory system is about 82% at 160 GeV/c. This polarisation is achieved 
in an elegant way, resulting from parity violation decays of the mesons ir + and K + , 
which in turn are produced by collisions of the SPS proton beam on a thick absorber. 
A sketch of the M2 beam line is displayed in Fig. 1. 



Naturally polarised muon beam: P (l = 82% 




Figure 1: Sketch of the production of a polarised muon beam for COMPASS. Note 
that the muons are 100% polarised in the pion rest frame (due to helicity conserva- 
tion). These polarised leptons, fi + , are the probes used to scan the gluon polarisation. 

The COMPASS target used in 2002, 2003 and 2004 consists of two cells contain- 
ing polarisable material in the solid state, each 60 cm long and 3 cm in diameter, 
separated by a gap of 10 cm. For the study of the gluon polarisation, the cells are 
longitudinally polarised in opposite directions (in order to allow for the extraction of 
spin asymmetries). The experimental acceptance for the resulting events is limited 
to an angular range of ±70 mrad. However, in the years of 2006 and 2007 the target 
setup was upgraded to a 3-cell system with an angular acceptance of ±180 mrad. 
These two improvements represent the only differences compared to the target setup 
used until 2004. In 2006 and 2007 the upstream and downstream cells are 30 cm long, 
and are both 5 cm away from a central cell of 60 cm long (their diameter is 4 cm). 
The former cells are longitudinally polarised in the opposite direction than that for 
the middle cell. 

Solid state targets are used, for all years of data taking, in order to achieve the 
high luminosity required for the physics programme of COMPASS. The target mate- 
rial used to study the spin structure of deuterons and protons is 6 LiD (2002-2006) and 
NH 3 (2007), respectively. The target cells are surrounded by a powerful supercon- 
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ducting solenoid which is used for their polarisation. The solenoid aperture defines 
the acceptance for event detection in the COMPASS spectrometer. The average po- 
larisation for the 6 LiD material is P t ~ 50%, with a dilution factor / 0.4. The 
latter accounts for the fraction of polarisable material inside each target cell. For the 
NH3 material, the average polarisation is Pt ~ 90% and the dilution factor is / ~ 0.15. 

The particles resulting from inelastic scattering collisions, i.e., the 160 GeV/c muon 
beam scattered off the 6 LiD or NH 3 material of the COMPASS target, are then 
detected in a two stage spectrometer whose detailed description can be found in 
Ref. [16]. Basically, each spectrometer is defined by a large magnet surrounded 
by several detectors (tracking, calorimeters, etc.). The first one is the large angle 
spectrometer, dedicated to the reconstruction of low momentum tracks, whereas the 
second spectrometer is designed to detect particles at smaller angles. A sketch of the 
COMPASS spectrometer is displayed in Fig. 2. 



|x filter . 




Figure 2: Artistic view of the COMPASS spectrometer during the 2004 data taking. 
The polarised muon beam enters from the left, and after being tracked by a set of SciFi 
and Silicon detectors collides in the polarised target. The products of the reaction 
are detected by a two stage spectrometer, using the SMI and SM2 magnets. 



The most common charged particles (e, tt, k and p) can be identified by a Ring Imag- 
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ing CHerenkov (RICH) detector, which is located in the first spectrometer. It consists 
of several multi-wire proportional chambers, containing Csl photocathodes, which de- 
tect the UV components of the Cherenkov light. In 2006, there was a considerable 
upgrade of the RICH: the central part of the detector was replaced by multi-anode 
photomultiplier tubes, yielding a larger number of photons detected (by extending 
the range of detection to the visible wavelength) together with a much faster response 
(reducing the uncorrelated background of halo muons). Together with a readout elec- 
tronics refurbishment, this upgrade allowed for a much cleaner particle identification. 



3 Event selection 

The cleanest way to measure the gluon helicity in the nucleon is to analyse the events 
from polarised inelastic scattering experiments. In order to be sensitive to Ag(x), 
one must tag a process involving a polarised lepton-gluon interaction. In COMPASS, 
one of the possibilities to do it is to reconstruct charmed mesons: 




Figure 3: Open charm production from a polarised PGF process (at LO) 



This method provides a clean signature of PGF events, because the intrinsic charm 
content of the nucleon is negligible in the kinematic domain of COMPASS (cf. Fig. 
4). In this analysis the PGF process is tagged using the D° meson production. This 
production is limited to a range of small XBj, XBj < 0.1, because at COMPASS kine- 
matics the cross-section for PGF events decreases rapidly with xbj- 
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Figure 4: Intrinsic charm predictions obtained for the virtual photon energy v equal 
to its average value for the COMPASS data: (u) = 0.95 GeV. The experimental 
points represent measurements of the charm structure function by the collaboration 
EMC [17], and the black dotted curves are theoretical fits for the PGF process using 
different renormalisation and factorisation schemes [18]. The blue curve shows the 
intrinsic charm prediction assuming a nucleon fluctuation into charmed meson-hadron 
pairs [19]. The red dashed curve shows the intrinsic charm prediction assuming charm 
fluctuations at the partonic level [20]. 

The COMPASS spectrometer is designed to reconstruct the D° mesons through the 
invariant mass of their decay products, Kit pairs, and for that purpose the RICH 
detector plays an extremely important role: the requirement of proper identification 
of kaon and pion candidates reduces significantly the combinatorial background un- 
derlying the resonance, which is centred on the D° mass, cf. Fig. 5. In addition, 
the following kinematic cuts are applied: on the fraction of the virtual photon energy 
carried by the D°, 0.2 < Zpo < 0.85, and on the angle between the charmed meson 
direction and the resulting kaon in the D° centre-of-mass system, |cos#*| < 0.65. 
They are important to reduce the contamination of the PGF sample by events com- 
ing from processes that involve the fragmentation of a struck quark, because most 
of these events are almost collinear with the virtual photon direction or have a zjya 
values close to zero. The combinatorial background can be further reduced using the 
following channel: D* — > D°n s i ow with D° — > Kit (D° tagged with a D*). By apply- 
ing a cut on the difference of the reconstructed masses for the D* and D° mesons, 
3.2 MeV/c 2 < M%? - M T D f — M w < 8.9 MeV/c 2 , one can verify that there is not 
much room left for the slow pion momentum. Since the above mass difference is re- 
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constructed accurately 2 , it can be used as a cut to significantly improve the purity of 
the D° signal. Consequently, three new channels of lower purity can also be studied 
(all tagged with a D*): D° — > K7nr°, D° — > K su t,iT, and D° — > Kmnr. The resonance 
observed in the first channel emerges from the combinatorial background as a 'bump', 
centred around -250 MeV/c 2 in the D° — > Kn spectra, due to the fact that the extra 
7T° is not directly reconstructed for this analysis. The second channel represents the 
D° candidates without RICH response for the kaon mass hypothesis (sub-threshold 
kaons with p(K) < 9 GeV/c), and the last one involves another D° decay mode 
which helps to improve the statistical precision of the Ag/g measurement. The final 
D° samples used in this analysis are shown in Fig. 5. 



4 Method to extract the gluon polarisation 

The number of D° candidates collected in a given target cell and time interval is: 




A^' D ° X + -\a b 
b s + b 



(7) 



^n^h'd"x j g longitudinal double spin asymmetry of the differential cross-section 
for events with a D° or D in the final state, and Ab is the corresponding asymmetry 
originating from combinatorial background events. Furthermore, m = M^n (or m = 
MxTrnn) and the symbol £ denotes a set of k — 1 kinematic variables describing an event 
(Pt° , E D o, ...), whereas a, and i] are the spectrometer acceptance, the incident muon 
flux integrated over the time interval, and the number of target nucleons respectively. 
The differential unpolarised cross-sections for signal and background events folded 
with the experimental resolution as a function of m and £ are represented by s = 
s(m,£) and b = 6(m,£) respectively. The ratio s/(s + b) represents the signal purity 
of the D° mass spectra. The information on the gluon polarisation is contained in the 
muon-nucleon asymmetry A^ N . The latter is defined in LO QCD as a convolution 
between the ratio of polarised/unpolarised partonic cross-sections, (Aa/a)^ 9 ^ 11 cc , 
and the ratio of polarised/unpolarised gluon structure functions (Ag/g): 



A» N = (a LL }^_ with a LL =^={ ^ & k \. (8) 



Four equations like Eq. (7) are defined, i.e. one equation for each cell and spin 
configuration of the target. Since the factors s/(s + b) and a LL (= ((III)) have a large 
dispersion, a weighting method is used to minimise the statistical error. The signal 
weight, ojg = P^fdLL[ s /( s + b)], and the background weight, ujb = P^fD\b/(s + b)], 



2 Note that the cryogenic solid state target of COMPASS does not permit the distinction between 
the D* and D° decay vertices (or between the primary and the decay vertices of a 
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Figure 5: Invariant mass spectra for the D* Knr , D^- n7T o, D* KlTV1T , D* K and D° Kv sam- 
ples. The symbol D* denotes a D° meson tagged with a D*. The subscript notation 
indicates the final state of a D° decay. Both deuteron and proton data are included 
in these samples. 
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are used, where D represents the transfer of polarisation from the muon to the virtual 
photon. By weighting the events from each of four equations (like Eq. (7)) with u s , 
and then similarly with wg, it becomes possible to solve the system of 8 equations in 
order to extract Ag/g 3 . Only 7 unknowns exist in this system of 8 equations: Ag/g, 
As and 5 independent quantities involving the product acpr] (see Ref. [21]). How- 
ever, to solve the system in an optimal way one needs an and s/(s+b) for every event. 

In the next-to-leading order (NLO) approximation, the physical asymmetry is decom- 
posed as follows: 

A*» = a% F ^ + al L A 1 , (9) 

where in a^ F all the soft, virtual and gluon bremsstrahlung corrections to the LO- 
PGF process are included (cf. diagrams a), b) and c) of Fig. 6). The asymmetry 
a\ L accounts for the presence of physical background inside the D° resonance. It 
corresponds to an interaction between the virtual photon and a light quark, which 
subsequently (or immediately before) emits a gluon that produces a cc pair (cf. dia- 
gram d) of Fig. 6) . One of these charmed quarks may fragment into a non-PGF D°. 
Finally, A 1 represents the inclusive asymmetry which is experimentally well known 
[22]. The contamination described by the second term of Eq. (9) is negligible and it 
appears only at NLO or higher orders of approximation in QCD. 

Using Eqs. (7) and (8) or Eqs. (7) and (9) one is able to determine Ag/g at the 
LO or NLO approximations. In order to make possible the use of COMPASS results 
in global analyses, a set of virtual photon asymmetries, A 1 * N ~*' IJ,D x = A tlN /D, was 
determined in bins of and Ejjo. They are obtained from Eq. (7) using as a signal 
weight the factor u s = P^fD[s/ (s + b)]. Note that determined values of A 1 * 1 ^ are in- 
dependent of ail and, therefore, their determination does not depend on a theoretical 
interpretation. The criterion for the choice of the binning is to provide asymmetries 
which are independent of the COMPASS acceptance for the PGF process. A table 
containing values of A^* N (p^° , E D o) can be found in Ref. [21]. 



5 Analysing power 

The analysing power for the open charm production, a^, is dependent on the full 
knowledge of partonic kinematics. Consequently, this asymmetry is not experimen- 
tally accessible because only one D° meson is reconstructed per event; the information 

3 COMPASS extracts Ag/g instead of Ag because experimentally it is much simpler to measure 
asymmetries than cross-section differences. 
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(a) (b) 




Figure 6: Examples of the NLO processes contributing to the muoproduction of the 
cc pair: a) virtual correction, b), c) gluon bremsstrahlung, d) light quark background. 
The remaining diagrams can be found in Ref. [23]. 

associated with the second charm quark is lost. Nevertheless, o,ll can be obtained 
from a Monte Carlo generator for production of heavy flavours (AROMA) . 

5.1 The analysing power at LO in QCD 

The AROMA generator is used without parton showers in order to generate D° 
events from LO-PGF processes. After a full Monte Carlo chain, where the generated 
events are constrained to the COMPASS acceptance, all D° mesons are reconstructed 
in the simulated spectrometer. This model reproduces very well the measured D° 
kinematics at COMPASS [24]. The polarised and unpolarised partonic cross-sections 
are calculated using the information from the generator, and then the kinematic 
dependences of aiL are parametrised with the help of a Neural Network [25] . Finally, 
the Monte Carlo parametrisation is used to obtain an for each real data event. 

5.2 The analysing power at NLO in QCD 

The phase space needed for NLO real gluon emission processes, 7*g — > ccg, is simu- 
lated through parton showers included in the AROMA generator. For every simulated 
event, an integration is performed over the energy of the unobserved gluon in the NLO 
emission process. The details about the upper limit of integration can be found in 
Ref. [21]. This integration reduces the differential cross-section for a three-body final 
state (ccg) to a two-body one (cc), which can be combined with the LO cross-section 
(cc, PGF) and the two-body virtual and soft NLO corrections. The procedure guar- 
antees a correct infra-red divergence cancellation. In this way, the total partonic 
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cross-section at NLO is calculated on an event-by-event basis for the spin averaged 
as well as spin dependent case, and consequently a^ F at NLO accuracy is obtained. 
The same procedure is applied for the correction originating from a q LL . In the NLO 
approximation, Ag/g is estimated from A"'* N using (a^ F ) and (a 9 LL ). In Fig. 7 one 
can see that the NLO corrections are important for this open charm analysis. 




Figure 7: Distribution of the analysing power an at LO and at NLO approximations 
for the simulated D* K ^ events. 



6 The signal purity 

The Neural Network (NN) described in Ref. [25] is also used to parametrise s/(s + b) 
on real data. Here, the goal is to obtain D° probabilities for every event. In Fig. 
8 one can see the outcome of this parametrisation: the D° mass spectrum reveals 
a probability behaviour in bins of [s/(s + &)]nN' i- e - its purity increases towards 
[s/(s + &)]nn = 1" Consequently, the statistical precision of Ag/g is significantly 
improved due to a good separation of the physical events from the combinatorial 
background. To accomplish this, two data sets are used as inputs to the Network. 
The first one contains the D° signal and the combinatorial background events. These 
events are called good charge combination ones (gcc), referring to the charges of par- 
ticles from D° decays, and they are selected as described in Sec. 3. The second set, 
the wrong charge combination events (wcc), is selected in a similar way except that 
the sum of charges of corresponding particles should be different from zero. It contains 
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Figure 8: The Kn invariant mass spectra in bins of the NN signal purity, [s/ (s+&)]nn> 
for the D* Kn sample (example from 2006 data). The last panel shows a comparison of 
the two purities, ([s/(s + &)]nn) an< ^ i s /( s + ^)lfit' ^ e ^ a ^ er ^ s obtained from a fit 
to the mass spectra. Curves show the background component of the invariant mass 
fits. The significance of the D* K7T signal is shown as the ratio s/b. 

only background events and is used as a background model. The NN performs a multi- 
dimensional comparison of gcc and wcc events in a ±40 MeV/c 2 mass window around 
the D° mass. Within the gcc set, signal events are distinguished from combinatorial 
background by exploiting differences between the gcc and wcc sets in the shapes of 
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distributions of kinematic variables as well as multi-dimensional correlations between 
them. An example of a properly chosen variable for the Network is the kaon angular 
distribution in the D° centre-of-mass system, as shown in Fig. 9. The distributions 
in the side band bins, shown in the right plot of Fig. 9, illustrate the good quality of 
the background model. 




Figure 9: Example of the gcc and voce distributions of |cos 6*\ in the D° centre-of- 
mass. Left / Right: region of the D° signal / side bands. 

The purity [s/( s + ^)1nN * s obtained from a simple function applied to the NN output 
(see Ref. [21]). Thereafter, the mass dependence is added as a correction from a fit to 
the D° mass spectra in bins of [s/( s + ^)]nN - By respecting the correct D° kinematic 
dependences, this parametrisation allows us to use s/(s + b) inside us in an unbiased 
way. 

7 Results and Conclusions 

The experimental results on the gluon polarisation are: 



/A \ L0 

(—) = -0.08 ± 0.21 (stat) ±0.09(syst) ®(x g ) = O.lllg;^. (10) 
/A \ NL0 

(—) = -0.20 ±0.21 (stat) ±0.09(syst) @(x g ) = 0.28lg;^. (11) 

Both results are obtained at a scale of (fi 2 ) = 13 (GeV/c) 2 . The LO value is compat- 
ible with all the world measurements of the gluon polarisation, as shown in the left 
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plot of Fig. 10. From this figure one can conclude that the gluon polarisation is small 
and compatible with zero within the range 0.07 < x < 0.20. Since the method used 
for the extraction of Ag/g at NLO is currently being improved, the result shown in 
Eq. (11) is still preliminary (the final result might be slightly different). Nevertheless, 
the presented value corresponds to the first world measurement of Ag/g at NLO. 

In the right plot of Fig. 10, the two open charm results for xAg are compared to 
the global fits: the experimental measurements favour small values of the gluon con- 
tribution to the nucleon spin, i.e. AG. Note that the positive solution 4 of COMPASS 
for AG (Jq 1 g(x)dx = 0.27 ± 0.09) is approximately 3cr away from the NLO point. In 
conclusion, the spin puzzle of the nucleon still remains to be solved. 
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Figure 10: Left: world results on Ag/g from LO analyses of the PGF process. Right: 
parametrisations of xAg(x, Q 2 ) together with the LO and NLO results obtained from 
the open charm analysis. The experimental points of the right plot are obtained from 
the well known unpolarised gluon distribution g, i.e. the Ag/g results are multiplied 
by g at the x points of the measurements. 
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